Large tunable photonic band gaps in nanostructured doped semiconductors 
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A plasmonic nanostructure conceived with periodic layers of a doped semiconductor and passive 
semiconductor is shown to generate spontaneously surface plasmon polaritons thanks to its periodic 
nature. The nanostructure is demonstrated to behave as an effective material modeled by a simple 
dielectric function of ionic-crystal type, and possesses a fully tunable photonic band gap, with widths 
exceeding 50%, in the region extending from mid-infra-red to Tera-Hertz. 
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I. INTRODUCTION 

The science of surface plasma waves on metallic- 
dielectric interfaces 1 , has received such interest as to 
become a true branch of physics, the palsmonics, and 
opened new perspectives in the control of light-matter in- 
teractions, see e.g. Refs^r— . The surface plasmon polari- 
ton (SPP), which results from coupling electromagnetic 
field with collective oscillations of electrons supported by 
the metal/dielectric interface, has unique physical prop- 
erties based on enhanced nanolocalized optical fields, al- 
lowed by the negative dielectric constant of the metal 
below the plasma frequency. For example SPP propaga- 
tion can be controlled by waveguides 5 or by plasmonic 
crystals whose metal film periodically nanostructured in- 
duces the plasmonic band gap 6,7 . 

Engineering of surface plasmons using nanostructura- 
tion made it possible to develop a new range of materi- 
als with remarkable optical properties such as extraor- 
dinary optical transmission 8 , optical filtering 9 , optical 
magnetism 10 , second harmonic generation 11 or higher 
harmonic generation in extreme UV— , and also to con- 
trol the optical processes at the femto-second scal o 13 i 14 . 
The main concern is to control the optical properties as- 
sociated with surface plasmons by engineering effective 
materials through periodicity, size and metal shape of 
the nano-objects. For instance a recent theoretical study 
proposed the nanosturcturation of metallic nano-sticks 
along the three direction of space to highlight a strong 
coupling between the incident light and the free electrons 
of the metal which opens a photonic band gap (with a 
stop band of 6.3%) in the telecom wavelength ranged. 

The main idea is to ensure opening of the gap by using 
the polaritonic nature of surface plasmons rather than 
Bragg reflections. As a matter of fact, the dispersion re- 
lation of an electromagnetic wave propagating along the 
interface between a metal and a dielectric, the so-called 
surface wave, is similar to the one resulting from a strong 
coupling between a photon and an oscillator, it is said to 
be ionic-crystal-likei. However, while the low energy sur- 
face plasmon branch is purely two-dimensional, the high 
energy branch of the dispersion relation is radiative and 
therefore a photon propagating along the interface will 
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not see the stop band. Studies of such systems dates 
back to 1969 16 , revealing three branches in the SPP dis- 
persion relation, two of them corresponding to the anti- 
symmetric modes, the third one to the symmetric modes. 

We study light propagation in a periodic nanostructure 
constituted of planar layers of a metal and a dielectric, 
sketched in FigJTJ assumed infinitely periodic and infinite 
in the second transverse direction. We demonstrate that 
this metal-dielectric nanostructure, worked as a planar 
waveguide array, but outside the guided mode regime, 
naturally and spontaneously generates SPP and then be- 
haves as an effective three-dimensional material whose 
optical properties are completely and precisely under- 
stood by means of a single effective dielectric function 
of the ionic-crystal type, the formula ([2]) below. 

Moreover, using, instead of a metal, a doped semicon- 
ductor is a means, by adapting the doping level, to set 
the plasma frequency in a chosen range, such as, in par- 
ticular, to avoid absorption due to interband transitions. 
Then the simplicity of the obtained effective dielectric 
function allows one to determine readily how the choice 
of the materials and the dimensions of the nanostructure 
will tune the width of the stop band, and how the doping 
level will tune its position. The properties of the effective 
material in the gap region is the principal concern of our 
study with future interest in its switching properties in 
nonlinear regimes, see e.g. Refs j 17 i 18 . 

We thus demonstrate the polariton nature of the sys- 
tem when it is worked with transverse magnetic (TM) 
fields, in the long-wave approximation. Similar studies in 
literature focus mainly on the SPP properties of metal- 
dielectric structures, looking thus at the lower branch of 
the dispersion curve. Our approach predicts the behavior 
of the system on the entire spectrum, namely the lower 
polariton branch, the stop band and the upper polariton 
branch. Last we show that the approach applies as well 
for transverse electric (TE) fields. The structure works 
then as a simple metal for which we compute explicitly 
the high frequency dielectric constant and the (tunable) 
plasma effective frequency. 
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II. MAIN RESULT 



From now on, all frequencies are normalized to the 
plasma frequency uo p the wave numbers to k p = uj p /c, 
the lengths to k~ x (including spatial variables) and time 
to ujp 1 . Moreover we name frequency the actual angular 
frequency, units and all dielectric constants are 

the relative ones. The doped semiconductor working as 
a metal, we use a Drude model with optical index n(uu) 
and thus the final index profile n(z,uj) of the periodic 
array can be written in the elementary cell as 



n 2 (z, oj) 



z€ [-6,0], 
z e [0,a]. 



(1) 



The damping coefficient 7 is also dimensionless and 



scaled to uj p . 
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FIG. 1. (color online) LEFT: scheme of the nanostructure, 
infinitely periodic in direction z and infinite in direction y. 
RIGHT: example of the SPP dispersion law oo(q). The pa- 
rameter values are given in Eq.JJ]). The red and blue curves 
are the true dispersion relation solutions of Eqs. (|1Q|) and 
respectively, the black curves (almost indistinguishable 
from the red curves) are the plot of the ionic crystal disper- 
sion law (|12|) . The dashed line is the light line in the passive 
semiconductor. 

We shall demonstrate that, under TM irradiation, in 
the long wave limit, the resulting metamaterial is equiv- 
alent to a single layer having the following effective di- 
electric function 



uj(uj + 27) — CJ 2 ' 



(2) 



with the resonant dimensionless eigenfrequency uj r and 
high frequency dielectric constant e given by 



be 



as, 



be 1 



(a + b)e 1 e 
~be~' 



as, 



(3) 



Eq.Q is identical to the dielectric function of an ionic 
crystal with relative high frequency dielectric constant e, 
characteristic transverse frequency uo ri and longitudinal 
frequency uo p = 1 (by normalization). As a matter of fact, 



Fig H] displays an example of the exact dispersion relation 
for the anti-symmetric modes (red curves) compared to 
the ionic crystal dispersion law uj = q/y/s eff (plotted with 
7 = as black curves). The blue curve corresponds to 
the symmetric modes as explained in SeclIVl 



III. SIMULATIONS 

We perform numerical simulations of a realistic sys- 
tem submitted to TM normal incidence and calculate the 
reflection and transmission coefficients. The presented 
spectra are simulated with the commercial software pack- 
age running finite difference time domain—. We shall 
illustrate our task by using InAs for both semiconduc- 
tors (high-frequency dielectric constants e 1 = e = 11.7), 
doping with silicon (Si) at a density of 10 20 that pro- 
duces a plasma frequency uo p = 3.42 10 14 rad/s and hence 



1.14 r ad I /ami*. The dimensions of the nanostruc- 



ture and the damping factor are taken as 

a = b = 0.2 fim, h = 1 /im, 7 = 10 13 rad/s, (4) 



which means normalized quantities a 
1.14 and 7 = 0.03. 



0.23, h 
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FIG. 2. (color online) Plots of the reflectance (black curve) 
and transmittance (red curve) obtained (a) from simulations 
of the structure at normal incidence, and (b) by the effec- 
tive dielectric function J2} with formulas f|13|) . The resonant 
frequency uj r given in Q defines the stop band [cj r , 1]. 
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Fig|2] shows the result of a simulation where transmit- 
tance and reflectance are plotted in terms of the fre- 
quency of the incident TM plane wave at normal inci- 
dence. It is then quite remarkable that the spectra are 
almost identical, which evidences our main claim, though 
Figj2jD has been obtained with two major approxima- 
tions: the ionic-crystal dielectric function (J2j) and the 
effective layer modelization. 

The obtained normalized stop band [uo ri 1] can then be 
easily broadened by playing with the dimensions of the 
nanostructure, the plasma frequency uo pi and the ma- 
terials. For example changing the passive semiconduc- 
tor from InAs to GaSb (dielectric constant e 1 = 14.4) 
increases the stop band from 29% to 33% for a depth 
h = 1 /im. Last, in order to confirm that the stop band is 
not the result of a Bragg scattering, we made a series of 
simulations with an incident angle varying from to 15 
degrees, and with layer thicknesses up to a = b = 0.7 /im, 
without any noticeable change in the spectra. 
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FIG. 3. (color online) Reflectance (black curve) and transmit- 
tance (red curve) obtained as in Fig(2]now with a = 0.3 \im 
and b = 0.1 /ira, the other parameters being the same. 

Another example is displayed in Fig|3] where we have 
modified only the geometry of the structure by increas- 
ing the width a of the doped semiconductor to 0.3 /im 
and reducing the width b of the passive semiconductor to 
0.1 /im. The size of the stop band then reaches 50% of 
uj p . Small discrepancies appear in the lower end of the 
spectrum, related to the long- wave approximation which, 
in the cases a ^ 6, is less accurate. We found that the 



lower branch of the dispersion relation (fT2j) is for a ^ b 
less close to the actual branch given by the solution of 

sun- 

It is now worth understanding the field intensity dis- 
tribution for some representative points of the spectrum. 
In the lower branch of the dispersion curve, the graph 
(1) of Fig|4] shows that transmission is accomplished by 
means of surface plasmons (in that case both k and k\ are 
pure imaginary numbers) and that the peak of transmis- 
sion is due to a stationary mode along Ox pinned at the 
interfaces and living essentially inside the passive semi- 
conductor. The next graph (2) is plotted for a frequency 
inside the stop band and shows a field exponentially de- 
creasing along Ox, as indeed q is now a pure imaginary 
number. The last graph (3) is taken for a frequency in 
the upper branch of the dispersion curve. It shows a 
transmission essentially inside the doped semiconductor 
(again with a stationary mode along Ox). This transmis- 
sion can be understood by noting that, at this frequency, 
the index of the doped semiconductor is smaller than the 
index e 1 of the passive semiconductor. 
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FIG. 4. (color online) Field intensity plots for the parame- 
ter values used in Fig(T] reproduced in the first graph where 
numbers indicate the 3 values of the frequencies which then 
index the corresponding intensity plots. 

With all these simulations in hand, it is time to demon- 
strate now the essential result, namely to derive the ex- 
pression of the effective dielectric function ((2]). This is 
done in the following section simply by using Maxwell's 
equations and continuity conditions for an infinitely pe- 
riodic structure. 



IV. TM FIELD: THE POLARITON MACHINE 

We thus calculate the dispersion relations and the re- 
sulting reflection coefficient by assuming the nanostruc- 
ture to act as a single layer with an effective dielectric 
function. This is done to show first that the optical 
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properties of the metamaterial indeed result from SPP 
generation. Second, it allows us to study the long wave 
limit and demonstrate that the dielectric function has the 
explicit simple approximate expression ((2]). 

Considering an infinitely periodic nanostructure, we 
may extract the elementary cell z G [— 6, a] and connect 
the interfaces z = — b to z = a by usual continuity rela- 
tions. As we expect to excite surface waves at the inter- 
faces, we first consider TM field and seek solutions 



E(x, z,t) 



E(z) 


F(z) 



J(ujt-qx) 



(5) 



as the continuity of the tangential components readily im- 
plies that the x-component q of the wave vector assumes 
the same values in either region. In the Drude model, the 
Maxwell equation reduces to the usual Helmoltz equation 
for E(z), together with an explicit expression for F(z), 
with the optical index function (pQ), namely 



d 2 E 

dz 2 



+ [n 2 (z,uj)uj 2 - q 2 ]E = 0, 
dE 



[n (z,uj)uj -q ]F = 



(6) 
(7) 



According to Eq. (pQ) we may define the eigenvalues k\ and 
k in each region by 



z G [-6,0] : 
z G [0,a] : 



7 2 2 2 

k 2 =eu 2 (l-^)-q 2 , 
ujuj' 



(8) 



where we define uj' = uj + ij. Solutions to (j6j) are then 
sought under the form 



z G [-6,0] 
z G [0,a] 



Ex(z) = A ie iklZ + B ie - iklZ , 
E(z) = Ae ikz + Be~ ikz . (9) 



Continuity relations in z = and z = a, —b for E(z) 
and F(z) then result in an algebraic homogeneous lin- 
ear system whose solvability condition produces a 4 x 4 
determinant. It actually factorizes in two terms which 
can be directly retrieved by assuming first the sub-case 
E(a) = —E(0), which gives 

e,k tan(afe/2) + eki(l — ) tan(6fci/2) = 0. (10) 

UJUJ 

This case is called anti-symmetric because, due to re- 
lation JJJ), the transverse field component F(z) is anti- 
symmetric on the elementary cell [—6, a]. Then the sym- 
metric case is obtained for E(a) — E(0) for which 

skM tan(afc/2) + e,k tan(6fci/2) = 0. (11) 

ujuj' 

Note that the complete electric field is neither symmet- 
ric nor anti-symmetric because of relation (0), still the 
dominant term is played by the transverse component 



F(z) which has a defined symmetry. The above two 
equations must be read with the definitions flS) as re- 
lations between uj and q. As usual these are implicit 
expressions solved numerically which, for the dispersion 
relation (j(q), furnishes the graph of Fig JT] where the red 
curve corresponds to the modes ([TO]) and the blue curve 
(almost a straight line) shows the modes obtained from 
Eq.JTTJ. These graphs are obtained as usual by assuming 
7 = to deal with real- valued expressions. 

We consider now the long-wave limit where a and 6, 
which are scaled to /c" 1 , are small quantities such as 
to allow the Taylor expansions tan(afc/2) ~ ak/2 and 
tan(6fci/2) ~ bk\/2. Note that such approximations are 
valid in a finite wave number domain. Then the disper- 
sion relation ([TO]) written for the wave number q is easily 
demonstrated to become, with help of definitions (j8]), 



„ 2 uj(uj^ij) 
euj 1 



uj(yj + ij) — uj 2 



(12) 



with the dimensionless resonant frequency uj r and dielec- 
tric constant e given in (|3]). The same treatment being 
applied to relation ([TT]) shows that the wave numbers 
k and ki cancel out and the relation simply reduces to 
ujuj' = uj 2 where uo t is defined below in Eq. ([2T|) . and which 
for a = b reduces to uj r . This explains why the blue curve 
of FigJT] tends to the straight line uj r (within the chosen 
wave number range). 

Therefore we can state that the nanostructure consid- 
ered here behaves, at least in the long- wave limit, as a sin- 
gle layer of active medium with effective dielectric func- 
tion £ eff , which from ([T2]) gives our fundamental result 
(|2]). This result can be used to evaluate reflectance and 
transmittance of a layer of depth h as 21 



R(uj) 
T(uj) 



7*1 + V2 ex.-p[2iujhn eii ] 



n 



T\T2 exjp[2iujhn ei 
t\ti exjp[iujhn e 



1 + rir2 exp[2iujhn e 



(13) 



with the following definitions (Fresnel formulas) at nor- 
mal incidence 



ri 



h = 



n 



i e{{ + n 
2no 



T2 



n3 



n 



n 



n 



t 2 = 



n 3 + n e 

2 ^eff 

n 3 + n ef 



(14) 



when the layer lies on a substrate with optical index 713, 
when no is the index of the medium in the incident re- 



gion and where n eff 



is given by the approximate 



expression (j2j). The graphs (b) of Figs. [2] and [3] show 
the plots of the above reflectance and transmittance ([T3]) 
with no = 77,3 = 1. Comparison with corresponding 
graphs (a) shows the accuracy of such a simple mod- 
elization. 

One may also compute from the relations ([TO]) and ([TT]) 
the dispersion laws uj(k) inside the doped semiconductor 
and uo(ki) outside. In the long- wave approximation, one 
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obtains that A: is a pure imaginary number outside the 
stop band [cj r ,l] while k\ is a pure imaginary number 
below the resonant frequency uo r . This is why the lower 
branch of the dispersion relation uj(q) < uj r represents a 
surface plasmon: the field is exponentially located at the 
interfaces. 



V. TE FIELD: A SIMPLE METAL 

It is interesting to consider the case of incident trans- 
verse electric (TE) fields for which we demonstrate now 
that the system behaves as an effective metal with a 
Drude-like dielectric function. Therefore the system il- 
luminated with TE fields does not give rise to polaritons, 
i.e. does not behave as a ionic-crystal, which is related 
to the fact that generation of surface waves on a metal- 
dielectric interface occurs only with TM fields. 

The Maxwell equation for the following TE field struc- 
ture (again the wave number q in the x-direction is the 
same in both regions by continuity) 

/ \ / \ 

Ei = Ei{z) e i M-^) ) e= E{z) ^M-^), 

V I V I 

(15) 

reduces to the Helmoltz equation (j6j). In that case the set 
of continuity relations has to be completed with the mag- 
netic field calculated by means of VxE = — <9tB, namely 




i(ut-qx) 



(16) 



and similarly for Bi. As in the preceding section, a solu- 
tion is sought under the expressions (|9]) but now the two 
sets of continuity relations at each interface are those ob- 
tained from (E — Ei) x n = and from (B — Bi) x n = 0, 
where n is the unit vector in the direction Oz normal to 
the interfaces. These sets are thus 

E(0) = E 1 (0), E(a)=E 1 (-b), 

d z E(0) = a,£i(0), d z E{a) = dzE^-b). (17) 

It is then a simple task to express the above relations and, 
as before, consider separately the two sub-cases E(a) = 
E(0) and E(a) = —E(0) for which we obtain respectively 

k tan(afc/2) + h tan(6fci/2) = 0, (18) 
ki tan(afc/2) + k tan(6fci/2) = 0. (19) 

The wave vectors k and k± are defined in Eq. (jSJ, there- 
fore the above relation furnish the dispersion relation 
u(q) in each case. 

We consider then the long-wave approximation for 
which the system is seen as a single layer, and for which 
the resulting effective dielectric function can be explic- 
itly computed. In the limit tan(afc/2) ~ ak/2 and 



tan(6fci/2) ~ bk\/2, the relation ([T8]) can eventually be 
expressed as 



2 2 

q = e n „(j , e. 



as + be 1 
a + b ' 



--^t 1 "^))' (20) 

(21) 



ae 



as + be 1 ' 



It is a Drude-like dispersion relation with a new and 
a new plasma frequency uo t (normalized to uj p ). Thus 
the effective equivalent material behaves just as a metal 
for TE modes. Reflectance and transmittance are given 
by formulas ([T3]) with, for normal incidence only, Fresnel 
expressions ([T4]) , and now with e eii = n? {{ given in Eq. 
([20]) . This result is compared with numerical simulations 
on Fig|5] with the same striking efficiency as in the TM 
cases. 
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FIG. 5. (color online) Reflectance (black curve) and trans- 
mittance (red curve) obtained with a = b = 0.2 \im for a 
TE incident field with parameters Q. The structure indeed 
works as a metallic effective layer with the effective dielectric 
function defined in Ea. ([20|) . 



Note finally that the long-wave limit for the second 
relation ([T9]) reduces to the two solutions k — and k\ = 
associated with the trivial dispersion laws obtained out 
of Eq. (j8j). Correspondingly, the TE field in that case 
results to be a constant (along z) and therefore simply 
vanishes due to the symmetry condition E(a) = —E(0). 
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VI. COMMENTS AND CONCLUSION 

When the long-wave limit is not accurate enough, one 
might think of solving numerically the dispersion rela- 
tions ([TO]) and (pT|) to calculate e eff . However, beyond 
the long- wave approximation, the main effect to take into 
account is the occurrence of guided modes, primarily in- 
side the undoped semiconductor. In that case one must 
reconsider the problem since the method of the effective 
layer would simply fail. 

The presence of two types of modes makes it uneasy 
to define in general the stop band which is not always 
[cj r ,l]. By looking at the dispersion relations, close to 
the long- wave limit, we found that the stop band is the 
interval [cj r , 1] when a < b, but becomes [cjt, 1] when 
a > b. However in that case, the interval [uj r ,uj t ] is filled 
with symmetric modes which are not excited under nor- 
mal incidence, the effective stop band being then still 
[o; r , 1]. Note that when a = 6, the 2 frequencies uj r and 
cut coincide, and the interval [uo r , 1] is a gap also at oblique 



incidence. 

By comparing realistic numerical simulations with the- 
oretical considerations, we have shown that the nanos- 
tructure described in Fig JT] works in TM regime as a sin- 
gle layer with the effective dielectric function (j2j), repre- 
sentative of SPP generation. This result provides a com- 
prehensive understanding of the system as a whole, on 
the entire spectrum, to show that it behaves as a ionic- 
crystal material. Moreover our result provides a powerful 
tool to explore the possibilities offered by such a meta- 
material, e.g. by varying the sizes, the materials and the 
doping level, with the opportunity to reach the visible 
range. 

Note finally that, in view of experimental realization, 
different techniques can be adopted as e.g., epitaxy, de- 
position, etching or implantation. Although these are 
well proved techniques, the resulting optical properties 
of the structure will be affected by imperfections such as 
interface roughness, diffusion, geometrical defects, which 
requires further studies. 
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